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The dwarf spheroidal galaxy Draco has long heen considered likely to be one of the brightest point sources 
of gamma-rays generated through dark matter annihilations. Recent studies of this object have found that it 
remains largely intact from tidal striping, and may be more massive than previously thought. In this article, we 
revisit Draco as a source of dark matter annihilation radiation, with these new observational constraints in mind. 

We discuss the prospects for the experiments MAGIC and GLAST to detect dark matter in Draco, as well as 
constraints from the observations of EGRET. We also discuss the possibility that the CACTUS experiment has 
already detected gamma-rays from Draco. We find that it is difficult to generate the flux reported by CACTUS 
without resorting to non-thermally produced WIMPs and/or a density spike in Draco’s dark matter distribution 
due to the presence of an intermediate mass black hole. We also find that for most annihilation modes, a positive 
detection of Draco by CACTUS would be inconsistent with the lack of events seen by EGRET. 

PACS numbers: 95.35.+d;95.30.Cq,98.52.Wz,95.55.Ka FERMILAB-PUB-05-538-A 

INTRODUCTION 

It has long been thought that dark matter particles could be observed indirectly by detecting the products of their annihilations. 
Such products, including gamma-rays, neutrinos and anti-matter, have been searched for using a wide range of experimental 
techniques jj]- Gamma-rays from dark matter annihilations, in particular, have been sought after using both satellite and ground 
based experiments. 

The potential astrophysical source of gamma-rays from dark matter annihilations which is most often studied is the central 
region of our galaxy. Recently, observations by the Atmospheric Cerenkov Telescopes (ACTs) HESS j3l, Whipple Q] and 
Cangaroo 0] have revealed the presence of a very bright gamma-ray source from this direction. The spectrum of this source 
has been measured in steadily increasing detail by the HESS collaboration 0. Although the first HESS data from this source 
was not inconsistent with a spectrum from annihilating dark matter 0, it is now becoming difficult to reconcile the HESS data 
with such a spectrum. Instead, it appears more likely that an astrophysical accelerator is responsible for this bright gamma-ray 
emission. As a result, future dark matter searches in this region will face a background that will be very challenging to overcome 
0 . 

Given this newly discovered background, it is important to consider other possible regions in which an observable rate of dark 
matter annihilation radiation may be generated. Such gamma-rays may appear as point sources external to our own galaxy, such 
as Andromeda (M31), M87 or the Large Magellanic Cloud 0, or as a diffuse spectrum generated by a large number of distant 
sources 00. Observable quantities of gamma-rays may also be generated in dark substructure within our own galactic halo. 
Again, this may appear as a diffuse spectrum from a large number of dark matter clumps , or may be dominated by a few of 
the most massive dwarf galaxies within the Milky Way, such as Draco, Sagittarius and Canis Major II llllTIl . 

In this article, we will discuss the prospects for detecting gamma-rays from dark matter annihilations in the dwarf galaxy 
Draco. We focus on this particular object for several reasons. First, of the most nearby and massive dwarf galaxies, the halo 
profile of Draco is the most tightly constrained by observations. Although other dwarfs may actually be brighter sources of 
dark matter annihilation radiation (this is likely for both Sagittarius and Canis Major fllll ). the rates from these objects cannot 
be estimated with as much confidence. Second, since dwarf galaxies are dark matter dominated, containing very few baryons, 
gamma-ray searches for dark matter in these regions are very unlikely to be complicated by the presence of astrophysical sources. 
In light of the challenges faced for dark matter searches in the galactic center, this is clearly an important consideration. 

A third reason that we chose to focus on Draco is the potentially exciting results of the CACTUS gamma-ray experiment. In 
recent conferences Ol, the CACTUS collaboration has stated that they have detected an excess of ^100 GeV gamma-rays from 
the direction of Draco. Although still preliminary, this result, if confirmed, would have dramatic implications for dark matter. 

The remainder of this article is organized as follows. In the following section, we calculate the annihilation rate of dark matter 
in Draco, and the resulting gamma-ray flux. We then discuss the prospects for MAGIC and GLAST to detect this flux, and then 
lastly turn our attention to the possible detection of Draco by CACTUS, and the implications of such an observation for dark 
matter. 
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GAMMA-RAYS FROM DARK MATTER ANNIHILATIONS IN DRACO 


Gamma-rays can be generated in dark matter annihilations through several processes. Most distinctive are those which result 
in mono-energetic spectral lines, XX 77^ XX or XX 7^- In most models, these processes only take place through 

loop diagrams, and thus the cross sections for such final states are quite suppressed, and lines are experimentally challenging to 
observe. 

A continuous spectrum of gamma-rays can also be produced through the fragmentation and cascades of most other annihila¬ 
tion products. The spectrum which results depends on the dominant annihilation modes. Parameterizations of the gamma-ray 
spectrum from dark matter annihilation can be found for several cases in Refs. QliQ. 

The normalization of this spectrum depends on the dark matter prohle of Draco. Assuming that the halo prohle is approximately 
spherically symmetric, the total annihilation rate in Draco within the radius Va is given by 
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where {aAv) is the WlMP’s annihilation cross section and is its mass. p(r) is the density of dark matter at a radius r from 
Draco’s center. This annihilation rate leads to an isotropic flux of gamma-rays from Draco that is given by 
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where D is the distance to Draco and is the number of gamma-rays produced per annihilation in the energy range of a give n 
detector. The distance to Draco has been determined to be 75.8 ± 0.7 ± 5.4 kpc from an analysis of RR Lyrae variable stars 11511 . 

In contrast to the halo prohle of the Milky Way galaxy, the properties of Draco’s prohle are somewhat constrained. For 
instance, it has been shown that Draco’s halo has not been tidally stripped due to the interaction with the Milky Way halo O. 
Despite these constraints, however, curi'ent observations have been unable to determine the slope of Draco’s inner halo prohle, 
being equally consistent with a ’’cusped” halo prohle, i.e. a density prohle ^ , with 7 between —0.5 and —1.5 as given by 

N-body simulations, or with a hat core (7 ^ 0). We will consider both of these possibilities. 

Considering an NFW prohle 
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where = rminirs and ya = ra/rg. 

Alternatively, we can consider a halo with a Hat central core, 
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In table 0 we give numerical values for the NFW and cored prohle cases, for some specihc values of j/a and ymin- These values 
allow us to rewrite Eqs. I^andl^ as 
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where A is the value for a given prohle found in table HI There are several things to note about these results. Firstly, the 
annihilation rate does not depend critically on the outer radius integrated out to. Varying the outer radius between the scale 
radius, 77 , and much larger values, the overall annihilation rate varies only by about a factor of 2 for a cored prohle, and much 
less for an NFW prohle. Secondly, unlike in the case of the galactic center halo prohle, the annihilation rate is not very much 
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Profile Type 

A{ra = rg) 

A{ra » Vg) 

NFW 

0.875 

1.0 

Core 

0.160 

0.323 

Cusp, 7 = 1.1 

1.29 

1.52 

Cusp, 7 = 1.2 

2.16 

2.63 

Cusp, 7 = 1.3 

4.03 

4.12 

Cusp, 7 = 1.4 

11.1 

12.5 

Cusp, 7 = 1.45 

25.7 

27.4 


TABLE I: Values of the parameter A, as used in EqlH for various halo halo profiles. In each case, rmin = 0 was used. See text for more 
details. 


lower for a cored profile than for the NFW case. In particular, the rate is reduced only by a factor of 3 to 5, is the same values 
of po and Vg are adopted. In addition to cored and NFW profiles, we have also shown in table |I]results for profiles with a denser 
cusp. In the case that the density in the inner halo scales with l/r~* rather than 1 /r as in the NFW case, we find somewhat larger 
annihilation rates, as expected. 

Up to these modest halo model-dependent variations in A, the annihilation rate and gamma-ray flux depends only on the 
quantity as well as the annihilation cross section and mass of the WIMP. Several observational constraints can be applied 
to Draco to constrain the quantities po Vg. Most important for our purposes are the constraints on the circular velocity, and 
the requirement that dark matter halos be abundant enough to account for the ~20 dwarf spheroidal galaxies in the local group. 
Also important are the constraints arrived at by the requirement that the first stars formed in Draco at least 10 billion years 
ago, and that Draco’s virial radius extends at least to the most distant stars (~1.2 kpc from Draco’s center). Collectively, these 
constraints limit Vg to values between 7 and 0.2 kpc, and po to values between 10^ and 10® M 0 /kpc® for the case of an NFW 
profile (the allowed ranges are somewhat smaller for the case of a cored profile) E. Interestingly, the allowed regions in the 
Ts-po plane are rather narrow strips, with a slope of roughly po oc Vg dl. This leads to little change in the quantity pgr® 
with variation of pg and Vg. For either an NFW or a cored halo profile, the allowed range of the quantity p^r® varies by only a 
factor of approximately 200 (0.03-6.3 x 10^® Mg/kpc®) M- 

Applying these constraints on pg and Vg as well as the allowed range of D, we arrive at the following maximal and minimal 
gamma-ray fluxes from Draco: 
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where Vg, is the number of gamma-ray produced per annihilation in the energy range of the detector. In each case an angular 
radius around Draco of 1° was considered (which impacts the value of A). 

This range of fluxes conservatively encompasses the range which can be expected from annihilating dark matter in Draco 
given the current constraints on its mass distribution. If less conservative profiles are considered, such as denser cusps, fluxes 
larger by 1-2 orders of magnitude are possible. A smaller annihilation rate and gamma-ray flux than the minimal cored profile 
case cannot be easily accommodated, however. ' 


^ Comparing our results to those of Evans, Ferrer and Evans El , their fluxes fall near the center of our (maximal to minimal) allowed range (in log units). 
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FIG. 1: The number of gamma-rays produced per dark matter annihilation for several annihilation modes. Results for annihilations to 
are shown as a dashed line, W^W~ and ZZ as dark and light (red) dotted lines, and bb and tt quarks as dark and light (red) solid lines. 


PROSPECTS FOR THE MAGIC AND GLAST EXPERIMENTS 


The experimental technology employed by the field of gamma-ray astronomy is currently developing very rapidly. This holds 
true for both ground-based Atmospheric Cerenkov Telescopes (ACTs) and satellite-based gamma-ray detectors. 

Of the currently operating ACTs, MAGIC is the best suited for observations of Draco. MAGIC’S northern hemisphere location 
at La Palma in the Canary Islands allows it to observe in the direction of Draco (unlike HESS, for example). VERITAS, once 
operational, will also benefit from its northern hemisphere location. MAGIC is also designed to have a lower energy threshold 
(~50 GeV for overhead sources) than other ACTs (~200 GeV for HESS or VERITAS), which is very important for dark matter 
searches. 

Erom the latitude of the MAGIC telescope (22° north), Draco reaches zenith angles as small as 29°. At this angle, an energy 
threshold of ~100 GeV should be possible. Large numbers of gamma-rays are produced above this energy only for WIMPs 
considerably more massive. In figure[2 the number of gamma-rays above 100 GeV per dark matter annihilation is shown as a 
function of the WIMP’s mass for several dominant annihilation modes. This clearly shows the difficulty in observing a WIMP 
not much more massive than the energy threshold of MAGIC, or other ACT. 

The primary background for ACTs is generated by hadronic cosmic rays. Eortunately, most of these showers can be identified 
and removed from the signal. In the energy range we are interested in, this background is roughly given by 




e X 10 ^ GeV ^ cm ^ s ^ sr ^ x 


/100GeV 

V Ebg 


2.7 


( 12 ) 


where e is the fraction of hadronic showers which are misidentified as electromagnetic, which is on the order of 1% for MAGIC. 
Integrating this above the 100 GeV threshold of MAGIC, and considering an effective area of ~ 5 x 10® cm^, this background 
accumulates at a rate of ~ 100 x e per hour over a 10“® sr solid angle (approximately a 0.1° by 0.1° circle). 

In the left frame of figure|2]we show the sensitivity of MAGIC to dark annihilations in Draco for three representative annihi¬ 
lation modes, and two halo profiles. We hnd that, for the case of a maximal NEW profile as discussed in the previous section, 
MAGIC will observe Draco with 5cr significance only for dark matter particles with annihilation cross sections of ~ 10“^® 
cm®/s or higher. This is somewhat larger than the maximum value for a thermal relic with a density equal to the measured cold 
dark matter density. Relics which are generated non-thermally may have larger annihilation cross sections, however. 

Eor the case of a satellite based experiment, such as GLAST, the prospects are quite different. GLAST is sensitive to gamma- 
rays down to 100 MeV, although to reduce the background, we will impose a threshold of 2 GeV, which is still far below that 
of MAGIC or other ACTs. With no background from misidentified hadronic cosmic rays, the diffuse gamma-ray background 
is all that needs to be overcome by GLAST. In the direction of Draco, this background has been measured by EGRET to 
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FIG. 2: The sensitivity of MAGIC (left) and GLAST (right) to dark matter annihilation radiation from Draco. The lower, darker set of lines 
adopts our maximal NFW profile, while the upper, and lighter, set of lines adopts our most conservative, minimal cored profile (see Has. 181111 . 
For MAGIC, we have considered a 5 x 10® cm^ effective area, a solid angle of 10“® sr, 40 hours of observation time, and 99% hadronic 
seperation (e = 0.01). For GLAST, we have considered a 10“* cm^ effective area, a solid angle of 9 x 10“® sr and one year of observation 
time. In each frame, the solid, dotted and dashed lines correspond to annihilations to bb, W'^W~ and T'^r~, respectively. All contours 
represent the cross section and mass required to generate a detection at the 5a level. 




be approximately 3.3 x 10“^ cm“^ sr“^ for gamma-rays above 2 GeV. Over one year of observation, a solid angle of 
0.3° X 0.3°, and a square meter effective area, this yields approximately 9 background events - a rate considerably lower that 
for an ACT. 

In the right frame of figure|2] we show the reach of GLAST to dark matter annihilations in Draco. GLAST clearly does much 
better than MAGIC for dark matter annihilating to heavy quarks or gauge bosons. In particular, thermally generated WIMPs 
with an annihilation cross section of ~ 3 x 10“^® cm^/s and lighter than ~500-700 GeV are detectable by GLAST in the case 
of a maximal NFW halo profile. The prospects are less promising for GLAST in the case of annihilations to taus, however, as 
few low energy gamma-rays are produced in tau decays. 

Some scenarios in which WIMPs are generated non-thermally can produce very large fluxes of gamma-rays from an object 
such as Draco. For example, in Anomaly Mediated Supersymmetry Breaking (AMSB) scenarios, the lightest supersymmetric 
particle is a nearly pure Wino which annihilates very efficiently to gauge boson pairs lIlTIl . If produced thermally, such a particle 
would only constitute a small fraction of the dark matter density. Non-thermal mechanisms can generate the observed relic 
density of Winos in AMSB scenarios, however Gi. 

In figure 13 we show the sensitivity of GLAST to neutralino dark matter in AMSB scenarios. Even in the case of the halo 
profile which produces the smallest allowed annihilation rate in Draco (the minimal core profile) GLAST will be very close to 
being capable of detecting gamma-rays from neutralinos in the ~80-200 GeV mass range. The absence of any signal from Draco 
seen by GLAST would thus exclude AMSB scenarios in this rather interesting mass range. If other halo profiles are considered, 
AMSB scenarios would likely be observable over a wide range of neutralino masses. 

In addition to the continuum signals discussed so far, dark matter particles can also directly generated gamma-rays through 
loop-diagrams, leading to mono-energetic spectral lines. If such a line could be detected, it would represent a ’’smoking gun” 
for dark matter annihilation. Such lines are very difficult to detect, however. For neutralinos, for example, the cross section for 
gamma-ray lines is no larger than « 10“^® cm®/s, and can be up to 4 to 5 orders of magnitude smaller Il4ll9ll . Also, only 1 
or 2 gamma-rays are produced in each of these processes. With these considerations in mind, we can consider an optimal case: 

= 30 GeV, {gav) = 10“^® cm®/s, = 2, and the maximal NFW profile. Together, this yields a flux of ~ 2 x 10“^^ cm“^ 
s“^, or approximately half an event per year of exposure for GLAST. While even one event at ~30 GeV would be intriguing 
from the direction of Draco, it would certainly not constitute a smoking gun of any kind. If a more steeply cusped halo profile 
were present, perhaps several line events could be seen by GLAST, but only in the most optimistic particle physics scenarios. 
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FIG. 3: The sensitivity of GLAST to neutralino dark matter in Anomaly Mediated Supersymmetry Breaking (AMSB) scenarios. The upper 
and lower dotted lines show the sensitivity of GLAST to dark matter annihilations for GLAST (for a 5a detection), adopting our minimal core 
and maximal NFW profiles, respectively (see Eos. 1811 iT . The solid line represents the cross section predicted for a nearly pure-wino dark 
matter particle, as is present in AMSB scenarios. Even with a halo profile which produces the minimum possible annihilation rate in Draco 
(the minimal core), winos with a mass of 80 to 200 GeV could be potentially detected by GLAST. With other halo profiles, much heavier dark 
matter particles can be discovered. 
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THE OBSERVATION OF THE CACTUS EXPERIMENT 

CACTUS is a ground based gamma-ray telescope located near Barstow, California. It is largely sensitive to gamma-rays 
above ^50 GeV, and has an effective area of up to ^50,000 m^ for ~TeV gamma-rays. CACTUS employs an array of mirrors 
which were designed for solar observations rather than for gamma-ray astronomy, and therefore is not as optimally suited for 
gamma-ray detection as other ACTs. It is hoped, however, that its enormous effective area will make up for these disadvantages. 

Very recently, the collaboration of the CACTUS experiment has, at a number of conferences, stated that they have detected an 
excess of gamma-rays from the direction of Draco 111 311 . Although this result is still of a preliminary nature, it is interesting to 
consider the implications of such a detection if it is confirmed. 

In 7 hours of observation, CACTUS detected approximately 30,000 events above background from Draco, of which roughly 
7000 and 4000 were above 100 and 125 GeV, respectively 11311 . To compare this with the predicted spectrum from annihi¬ 
lating dark matter, we convolve the injected spectrum with the energy dependent effective area of CACTUS, which has been 
parameterized as 11311 

Aeff «47,000m2[l-e-°°^^(®^-3® ®^®^) -F0.00025E.^(GeV)]. (13) 

Normalizing to the total rate seen by CACTUS, we can compare the energy distribution of events to that expected from annihi¬ 
lating dark matter. In table miwe show the number of events above 100 and 125 GeV expected for various dark matter masses 
and annihilation modes. The numbers reported by CACTUS appear to be consistent with the cases of a ^500 GeV dark matter 
particle annihilating to bb, a ~300 GeV dark matter particle annihilating to W~^W~, or a ~200 GeV dark matter particle anni¬ 
hilating to T+r”. We emphasize, however, that systematic uncertainties in CACTUS’S energy determination and understanding 
of backgrounds may modify these conclusions considerably. Despite these concerns, we conclude that annihilating dark matter 
with a mass in the 150-1000 GeV range appears to be consistent with the limited spectral information contained in the CACTUS 
signal from Draco. 

Each of these scenarios require very high annihilation rates in Draco, however, which leads us to two potential problems. 
Firstly, a very cusped or spiked halo distribution would be needed to accommodate this rate - roughly 10^ to 10^ times larger 
than the rate found for the maximal NFW model. To accommodate this, either a very large annihilation cross section (and a 
non-thermal production mechanism) or a very dense dark matter distribution (perhaps surrounding an intermediate mass black 
hole in Draco tifill ) would be required. 

Secondly, the EGRET satellite observed this region, and has placed limits on 1-10 GeV gamma-rays from this region of the 
sky. For most of the scenarios shown in the table, this limit is violated. 

Specifically, after imposing angular cuts of 1.71°, 1.18° and 0.82° for gamma-rays between 1-2 GeV, 2-4- GeV and 4-10 
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Total 

> 100 GeV 

> 125 GeV 

EGRET 

CACTUS Observation 

30,000 

7000 

4000 

- 

600 GeV, bb 

30,000 

9000 

5000 

290 

500 GeV, bb 

30,000 

7700 

3900 

400 

400 GeV, bb 

30,000 

6000 

2700 

630 

400 GeV, W+IV- 

30,000 

9200 

5100 

280 

300 GeV, W+IV- 

30,000 

7100 

3500 

470 

200 GeV, W+W- 

30,000 

4000 

1300 

1100 

300 GeV, r+r“ 

30,000 

15,000 

9500 

2.8 

200 GeV, r+r” 

30,000 

9200 

4200 

7.2 

150 GeV, r+r“ 

30,000 

5000 

1300 

16 


TABLE II: The approximate energy distribution of events reported by CACTUS compared to the prediction from various annihilating dark 
matter scenarios. The CACTUS observations appear to be consistent with a ~500 GeV dark matter particle annihilating to bb, a ~300 GeV 
dark matter particle annihilating to W^W~, or a ~200 GeV dark matter particle annihilating to t'^t~. In the last column, the number of 
events which EGRET should have seen is given for each case. 


GeV, respectively, EGRET actually observed 6 events, with an expected background of 4.1. Dark matter scenarios in the table 
which suggest that EGRET should have seen hundreds of gamma-rays from Draco are clearly inconsistent with this result. The 
exception to this problem, however, is for dark matter which annihilates to tau pairs. In this case, only a few events are expected 
to have been seen in EGRET. It is intriguing to note that this excess of ~2 events observed by EGRET, although certainly not 
statistically significant, is of the same order of magnitude as the rate expected for dark matter annihilating to r+r” given the 
spectrum reported by CACTUS. This annihilation mode may dominate, for example, in the case of a bino-like neutralino which 
annihilates through the exchange of a light stau. 

If the CACTUS signal is in fact a gamma-ray spectrum (as opposed to a poorly understood background, for example), GLAST 
will detect thousands of gamma-rays from Draco, determining its spectrum in detail. MAGIC should also easily detect such a 
source with high significance. With MAGIC currently operating, we expect that it will not be long before the CACTUS signal is 
either confirmed or invalidated. 


SUMMARY AND CONCLUSIONS 

In this article, we have revisited the possibility of detecting gamma-rays produced in dark matter annihilations in the dwarf 
galaxy Draco. Draco is the most well constrained of the Milky Way’s satellite galaxies, and therefore provides the best opportu¬ 
nity to make reliable predictions of dark matter annihilation rates and corresponding gamma-ray fluxes. 

Using the constraints on the dark matter distribution of Draco put forth in Ref. Ilbll . we have calculated maximal and minimal 
annihilation rates (and corresponding gamma-ray fluxes), considering both a cusped (NEW) profile and a profile with flat core. 
The variation that we And in the annihilation rate between even these two extreme scenarios is less than three orders of magnitude. 

We then proceeded to compare these rates to the sensitivity of MAGIC and GLAST. MAGIC is a currently operating ground 
based gamma-ray telescope, while GLAST is a satellite based gamma-ray detector scheduled to be deployed in 2007. We And 
that while both MAGIC and GLAST have the ability to detect dark matter in Draco in some scenarios {ie. a maximal NEW 
profile, low WIMP mass and favorable annihilation modes), dark matter in Draco can go undetected by these experiments in 
other cases. In some extreme cases (such as the non-thermal generation of neutralinos in anomaly mediated supersymmetry 
breaking models, for example), however, the lack of a detection by GLAST of gamma-rays from Draco could successfully rule 
out models, even if the most conservative halo model were assumed. 

Einally, we have also discussed the implications of the recent possible detection of Draco by the ground based gamma-ray 
detector, CACTUS. We And that to produce the signal reported by the CACTUS collaboration, an annihilation rate of dark 
matter in Draco is needed which is three to four orders of magnitude larger than can be accommodated for an NEW profile 
and an annihilation cross section consistent with thermally generated dark matter. Non-thermally produced dark matter and/or 
extremely high densities of dark matter in Draco would therefore be required to generate this signal. We also And that the 
CACTUS signal appears to be in conflict with the null results from the region by the EGRET experiment for most choices of 
the dark matter’s dominant annihilation modes. If the dark matter almost entirely annihilates to tau pairs, however, this conflict 
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can be (marginally) avoided. If CACTUS is in fact detecting this very large flux of gamma-rays from Draco, both MAGIC and 
GLAST should easily be able to confirm this result. 

We would like to thank Gianfranco Bertone, Joakim Edsjo, Francis Halzen, Peter Marleau, Karl Mannheim and Mani Tripathi 
for helpful discussions. DH is supported by the US Department of Energy and by NASA grant NAG5-10842. LB is supported 
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